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I.  Overview  of  the  Program 

The  propose  of  this  theoretical/computational  research  program  was  to  develop 
methods  and  to  investigate  the  fundamental  chemical  dynamics  of  elementary  processes 
imponant  in  systems  that  can  potentially  yield  large  amounts  of  energy  when  they  react. 

The  focus  was  on  fundamental  reactions  in  polyatomic  molecules  in  which  tunneling, 
electronic  nonadiabatic  processes,  and  vibrationally  mode-selected  reactions  occur  in  an 
effort  better  understand  some  of  the  basic  processes  and  some  of  the  specific  systems  of 
interest  in  developing  practical  high  energy  density  materials.  The  emphasis  of  the  research 
was  on  studies  of  molecules  with  four  or  more  atoms,  particularly  the  effects  of  many 
modes  and  the  energy  distribution  among  them  on  adiabatic  unimolecular  reactions, 
predissociation,  reaction  via  barrier  tunneling,  and  intramolecular  vibrational  redistribution 
(IVR). 

The  theoretical  methods  included  classical  trajectories  and  transition-state  theory. 
Semiclassical  corrections  were  employed  to  treat  nonclassical  processes  such  as  tunneling 
and  electronic  transitions.  We  have  extended  Monte  Carlo  transition-state  theory,i-5  which 
we  originally  used  in  treating  adiabatic  processes,  to  the  treatment  of  electronically 
nonadiabatic  (curve  crossing)  processes.^ 

II.  Publications,  Theses,  and  Oral  Presentations. 

The  following  is  a  list  of  the  publications  and  manuscripts  (in  preparation  for 
publication)  that  have  resulted  from  the  research  during  the  grant  period. 
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1918  (1992). 
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(1992). 
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193,347-354  (1992). 

Thomas  D.  Sewell.  Donald  L.  Thompson,  J.  Daniel  Gezelter,  and  William  H.  Miller,  "On  Some 
Problems  of  Correcting  the  Zero-Point  Energy  Problem  in  Classical  Trajectories  of 
Polyatomic  Molecules,"  Chem.  Phys.  Letters  193,  512-517  (1992). 

Karen  L.  Bintz,  Donald  L.  Thompson,  T.  R.  Gosnell,  and  P.  Jeffery  Hay,  "Mode  Localization 
and  Dissociation  Dynamics  in  the  CF2Br  Radical,"  J.  Chem.  Phys.,  in  press. 

David  Sahm  and  Donald  L.  Thompson,  "Semiclassical  Calculation  of  State  Specific 

Unimolecular  Electronic-Predissociation  Rate  Constants,"  Chem.  Phys.  Letters,  in 
press. 

Alison  J.  Marks  and  Donald  L.  Thompson,  "Classical  Trajectory  Study  of  Mode  Selectivity  in 
Deuterium-Subsrituted  Methyl  Isocyanide,"  in  preparation. 

Ward  H.  Thompson  and  Donald  L.  Thompson,  "A  Slater  Theory  Treatment  of  Electronically 
Nonadiabatic  Processes:  Transition  Rates  and  Mode  Mixing  Effects,"  in  preparation. 

David  Sahm  and  Donald  L.  Thompson, "  Unimolecular  Rate  Constants  for  the 

Predissociation  of  N2O:  Trajectory  Surface-Hopping  and  Phase-Space  Integral 
Evaluations,"  in  preparation. 

Yue  Qin  and  Donald  L.  Thompson,  "Tunneling  Effects  in  the  Cis-Trans  Isomerization  of 
Hydroxysilyene,"  in  preparation. 


The  following  theses  are  based  in  pan  on  work  supponed  by  this  grant: 

Thomas  D.  Sewell  (Ph.D.,  1991)  Thesis  Title:  Unimolecular  Reaction  Dynamics  of  Large 
Polyatomic  Molecules. 

Yue  Qin  (Ph.D.,  1992)  Thesis  Title:  Semiclassical  Tunneling  Corrections  for  Classical 
Dynamics  of  Unimolecular  Reactions  in  Polyatomic  Molecules. 

Karen  L.  Bintz  (Ph.D.)  Thesis  Title:  Theoretical  Investigations  of  Factors  Which  Affect 
the  Mode-Specific  Behavior  of  Unimolecular  Reaction  Dynamics. 


The  following  is  a  list  of  oral  presentations  of  the  research  results: 

"Mode  Specificity  in  Unimolecular  Reactions,"  Department  of  Chemistry,  Kansas  State 
University,  Manhattan,  Kansas.  May  24,  1990. 


") 


"Theoretical  Dynamics  Studies  ot  Unimolecular  Reactions."  Symposium  on  Compiiiaiional 

Chemisiry,  at  the  36ih  Annual  Pentasectional  Meeting  of  the  Oklahoma  Sections.  American 

Chemical  Society,  Lawton.  Oklahoma.  April  6,  1991. 

Lecture  Series  on  Chemical  Dynamics  at  the  Institute  ofTlieoretical  Chemistry,  Jilin  University. 

Changchun,  P.  R.  China,  August  16-23,  1991. 

III.  Review  of  Completed  Research 

The  research  focused  upon  the  funher  development  of  methods,  applications  of 
.some  existing  methods  to  new  problems,  and  studies  of  fundamental  dynamical  behavior  in 
systems  of  interest  to  the  KEDM  program. 

The  research  has  been  concerned  with  the  intramolecular  dynamics,  tunneling 
effects,  electronically  nonadiabatic  reactions  and  other  aspects  of  the  fundamental  processes 
that  affect  the  lifetimes  and  reactions  of  long-lived,  high  energy  systems.  The  central  focus 
of  this  program  are  the  dynamics  of  polyatomic  systems,  with  the  emphasis  on  realistic 
calculations. 

A.  Electronically  Nonadiabatic  Unimolecular  Reactions 

Nonadiabatic  processes  in  diatomic  molecules  are  well  understood  or  can  be 
straightforwardly  treated.  However,  curve  crossing  problems  in  many-mode  systems  are 
much  less  easily  treated  and  as  a  result  we  have  a  less  thorough  knowledge  of  the  dynamics 
of  them.  We  have  carried  out  several  calculations  to  study  the  dynamics  of 
predissociation^"'®  and  have  suggested  new  theoretical  approaches^dO  the  problem. 

We  have  used  the  classical  trajectory  surface-hopping  method*^  to  study  the 
nonadiabatic  spin-forbidden  dissociation 

N20(1I+)  N2(‘ V)  +  0(3P).  (Rl) 

We  chose  this  system  because  there  have  been  several  studies  of  other  aspects  of  it.  We 
used  the  model  potentials  of  Zahr  et  al.  The  crossing  seam  is  higher  for  bent 
configurations  than  for  linear  ones  in  this  potential.  We  investigated,  treating  all  degrees  of 
freedom,  the  effects  of  initial  mode  excitations  on  the  predissociation  rate.^  More 
specifically,  we  calculated  the  rates  for  excitations  of  the  bending  and  stretching  normal 
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modes.  We  also  studied  the  effect  of  the  choice  of  "hopping  coordinate"  and  the  influence 
of  the  weak-coupling  correction  on  the  rate  constant. 

The  rate  was  calculated  using  the  Landau-Zener  transition  probability  plus  the 
weak-coupling  correction's’ll  (when  the  Landau-Zener  probability  became  "large,  "  i.e.. 
greater  than  0.1).  About  10%  of  all  reactive  trajectories  involved  a  transition  requiring  the 
weak-coupling  approximation.  The  weak-coupling  correction  makes  a  significant 
difference  in  the  rates  (as  much  as  a  factor  of  2). 

.Although  it  is  understood  that  the  hopping-direction  coordinate  is  appropriately 
defined  as  the  direction  normal  to  the  seam,  we  also  calculated  rates  using  definitions  based 
on  internal  coordinates  (the  N2  center-of-mass  to  O  coordinate  and  the  shonest  N-O  bond  i 
since  they  are  easily  calculated  and  thus  are  sometimes  used.'^ 

One  of  the  main  interests  in  the  study  was  mode  selectivity  in  th.e  predissociation. 
The  rates  are  strongly  dependent  on  the  initial  excitation  because  of  incomplete  energy 
redistribution  along  with  a  significant  variation  in  crossing  point  energies  with  molecular 
geometry'.  That  is,  there  is  significant  changes  in  the  accessibility  of  the  "dissociation 
channel"  for  different  mode  excitations. 

The  calculations  just  described  were  done  using  "mode-specific"  initial  conditions. 
That  is,  normal  modes  were  initially  excited.  We  have  followed  up  these  calculations  with 
some  in  which  "state-specific"  initial  conditions  were  used.'-^  These  were  done  for  a  two¬ 
mode  collinear  model  of  N2O.  Assuming  that  most  trajectories  on  the  initial  potential  are 
quasiperiodic,  rate  constants  can  be  calculated  for  a  panicular  torus.  The  tori  were  chosen 
to  be  those  corresponding  to  semiclassical  states  according  to  EBK  theory.  As  would  be 
expected  on  the  basis  of  the  mode-specific  results,  there  are  wide  variations  in  the  state- 
selected  rate  constants.  These  calculations  confirm  our  initial  speculations  (made  when  we 
proposed  this  work  3  years  ago)  that  mode  specificity  in  curve  crossings  might  be  rather 
significant  due  to  the  required  "strong  directionality"  required  to  hop  from  one  surface  to 
another. 
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In  addition  to  these  dynamics  studies,  we  have  ak  developed  Monte  Carlo  phase- 
space  theory^  method  for  treating  nonadiabatic  unimolecular  reacrions.  In  this  approacli. 
the  rate  constant  is  nsidered  as  a  weighted  flux  through  a  "transition  state,"  which  is 
represented  by  the  locus  of  intersection  of  two  diabaiic  potential-energy  surfaces.  Phase 
space  points  in  this  region  contribute  to  the  rate  with  weights  given  by  the  transition 
probabilities  (we  used  the  Landau-Zener  equation  in  our  illustration  of  the  method).  The 
resulting  expressions  are  evaluated  by  using  Monte  Carlo  approximants.  This  leads  to  a 
general  method  that  is  easily  applied  to  large  systems  and  for  general  potentials.  We  have 
demonstrated  the  method  for  N2O  and  shown  that  it  gives  results  in  good  agreement  with 
classical  trajectory  surface-hopping  results.  Thermal  rate  constants  are  easily  computed 
using  this  method. 

An  important  pan  of  treating  nonadiabatic  reactions  is  the  accurate  representation  of  the 
potentials  in  the  region  of  the  seam.  Yarkonyl^  has  development  methc  obtaining  this 
information  from  ab  initio  quantum  mechanics.  Our  challenge  is  to  incorporate  this  information 
into  dynamics  calculations.  Usually  one  has  the  best  data  about  the  equilibrium  geometries  and 
(with  the  results  such  as  Yarkony  provides)  the  crossing  seam.  The  potentials  for  unimolecuUu- 
reactions  are  usually  fit  to  functions  based  on  valence  force  fields.  Thus  it  is  most  convenient 
to  write  the  potentials  in  the  crossing  regions  using  the  same  kinds  of  representations.  We 
have  explored  this  idea  by  using  Slater's  theory. 

One  of  the  strengths  of  Slater's  theory  of  unimolecular  reactions^®  is  that  it  provides  a 
way  of  visualizing  and  illustrating  some  fundamental  behavior  of  dissociating  molecules.  Its 
great  weakness,  of  course,  is  that  it  is  based  on  purely  harmonic  mode  motions.  Nevertheless. 
Slater's  theory  is  often  instructive  because  of  its  simplicity  and  because  it  relies  on  well- 
understood  notions  (i.e.,  normal  mode  motion)  to  treat  the  dynamics.  We  have  adapted 
Slater's  theory  to  treat  problems  involving  electronically  nonadiabatic  transitions. 

Slater's  approach  to  unimolecukir  reactions  involves  monitoring  the  value  of  a  time- 
dependent  critical  coordinate  qc,  which  is  expressed  in  terms  of  a  linear  combination  of  the 
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normal  mode  coordinates  Thus,  its  time  dependence  is  trivially  known.  Dissociation  is 
considered  to  have  occurred  when  qc  exceeds  some  maximum  value,  qo;  that  is,  when  there  is 
an  "up-zero  "  in  qc(t)  -  qo.  The  rate  is  then  related  to  the  frequency  of  up-zeroes. 

This  approach  has  been  generalized  by  Doll^*  and  applied  to  more  general  cases  of 
dissociation  as  well  as  surface  diffusion^^  and  desorption.^^  DqU  used  a  result  by  Kac^^  to 
obtain  an  expression  in  integral  form  for  the  frequency  of  up  zeroes.  Thus,  the  statistical 
propenies  of  qc(t)  can  be  utilized  rather  than  charting  the  coordinate  as  a  function  of  time  by 
calculating  classical  trajectories.  Doll  expressed  the  frequency  of  up  zeroes  in  terms  of  a  time 
average  and  showed  that  it  can  be  extended  to  an  average  over  phase  space.^' 

We  show  how  Slater's  theory  of  unimolecular  reactions  can  be  adapted  to  electronically 
nonadiabatic  transitions  and  provide  insight  into  some  fundamental  aspects  of  them.  The 
standard  semi-classical  approach  to  nonadiabatic  transitions  is  two-fold:  (1)  Classical 
Trajectory  Surface  Hopping  (CTSH);^^’^'*  and  (2)  statistical  methods.^'l^  in  the  first 
approach  a  classical  trajectory  is  evolved  on  the  initial  adiabatic  potential  surface.  Each  time  the 
trajectory  passes  the  crossing  seam,  defined  to  be  the  points  at  which  the  two  adiabatic  surfaces 
are  of  equal  potential,  it  is  given  a  certain  probability  of  "hopping"  to  the  other  surface,  which 
implies  completion  of  a  reaction  or  a  quenching  event.  The  hopping  probability  is  most 
commonly  determined  by  the  Landau-Zener  formula.  The  rate  is  calculated  from  an 
ensemble  of  trajectories.  This  is  costly  in  terms  of  computer  time  since  the  trajectories  in  the 
ensemble  must  be  followed  for  long  times.  Thus,  statistical  methods  are  attractive  since  they 
circumvent  this  problem. 

The  statistical  methods  generally  involve  averages  of  the  transition  probabilities  over 
phase  space.  Examples  of  such  approaches  are  given  by  Zahr,  Preston,  and  Miller,  Lorquet 
Cl  al.,—  and  Marks  and  Thompson.^  Zahr  et  al.  assume  that  the  trajectories  are  sufficiently 
long-lived  that  they  behave  ergodically.  This  assumption  allows  them  to  express  the  transition 
probability  as  an  average  over  phase  space.  These  methods  are  less  expensive  in  computer 
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time  than  the  CTSH  ap  oach.  However,  they  provide  much  less  insight  into  the  fundamental 
behavior  of  the  process. 

The  success  of  the  application  of  Slater’s  theory  to  the  problem  of  dissociation 
encourages  a  similar  treatment  of  nonadiabatic  transitions.  Anharmonicity  is  an  inherent  aspect 
of  dissociation  that  is  contradictory  to  the  assumption  in  Slater’s  theory  of  purely  harmonic 
behavior.  In  the  case  of  nonadiabatic  transitions,  one  is  interested  in  the  dynamics  which  take 
place  at  the  crossing  seam,  which  can  lie  relatively  low  in  the  potential  well.  As  pointed  out  by 
Lorquet  and  Leyh-Nihant,^^^  at  these  lower  energies,  anharmonicity  is  not  as  prevalent  and  the 
harmonic  approximation  more  closely  resembles  the  actual  potential.  Thus,  one  would  expect 
that  Slater  theory  might  be  well  suited  to  the  problem  of  nonadiabatic  transitions.  One  of  the 
primary  advantages  of  a  Slater  approach  to  nonadiabatic  transitions  is  that  it  provides  a  way  to 
obtain  information  about  the  mode  specificity  of  the  process.  Since  the  Slater  trajectory  is 
written  in  terms  of  the  normal  mode  coordinates  it  is  perfectly  suited  for  calculating  the  effect  of 
exciting  a  particular  mode  or  combination  of  modes  on  the  rate. 

In  summary,  we  have  showed  how  Slater  theory  can  be  adapted  to  treat  nonadiabatic 
transitions,  and  make  some  simple  applications  to  illustrate  it  usefulness.  Analytical 
expressions  for  the  crossing  seam  and  other  relevant  quantities  are  developed.  We  have 
described  how  Slater  theory  can  be  used  to  calculate  nonadiabatic  rates  for  real  systems  using 
ab  initio  information.  An  application  to  the  N2O  system,  using  a  model  potential,  that 
illustrates  the  method  (since  there  are  insufficient  ab  initio  results  to  fully  define  a  crossing 
.seam)  is  made.  We  have  also  shown  how  the  model  can  be  extended  to  include  anharmonicity. 
We  also  apply  Slater  theory  to  the  problem  of  mode  mixing  and  illustrate  how  it  can  model 
intramolecular/intermode  energy  transfer. 

B.  Tunneling  Corrections  for  Unimolecular  Reactions 

We  have  devoted  considerable  effort  to  developing  methods  and  carrying  out 
calculations  of  tunneling  in  unimolecular  reactions.^^'^^  Tunneling  may  be  an  imponant 
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factor  in  the  behavior  of  energetic  materials  if  they  involve  "storage"  in  a  hydrogen 
environment.  Thus,  it  is  important  that  we  be  able  to  treat  this  effect  if  we  wish  to 
accurately  model  the  behavior  of  these  materials. 

The  tunneling  splitting  of  the  eigenvalues  in  symmetric,  double-well  molecules  has  long 
been  a  topic  of  considerable  study.  The  symmetric  hydrogen-transfer  process  in  the 
ir'ilonaldehyde  molecule  (MA) 


has  been  extensively  studied.  Leopold  and  co-workers^®  reported  a  ground-state  tunneling 
splitting  value  of  21.58476  cm‘^.  Hutchinson^^  has  published  the  only  study  of  the  classical 
dynamics  of  hydrogen  transfer  in  MA.  He  employed  a  two  degrees-of-freedom  potential-energy 
surface  consisting  of  constrained  fhree-body  LEPS  and  htirmonic  angle-bending  terms.  Only  the 
motion  of  the  migrating  hydrogen  atom  was  treated  explicitly;  the  skeletal  C  and  O  atoms  were 
fixed  and  other  hydrogen  atoms  were  neglected  completely.  The  study  revealed  the  dynamics  to  be 
surprisingly  regular  even  at  energies  considerably  in  excess  of  the  barrier  to  hydrogen  transfer. 
Moreover,  the  classical  dynamics  are  relatively  insensitive  to  modest  variations  in  the  potential- 
energy  surface  parameters.  Tunneling  was  not  considered  in  the  calculation. 

We  have  used  classical  trajectories^^  in  conjunction  with  a  semiclassical  treatment  of 
tunneling33’34  to  calculate  the  splitting  in  the  Hutchinson  model.  We  used  the  classical-plus¬ 
tunneling  model  of  Waite  and  Miller^^  in  conjunction  with  a  more  recent  version  due  to  Makri  and 
Miller. Common  to  both  treatments  is  the  assumption  that  the  trajectory  evolves  classically  until 
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it  experiences  a  turning  point  in  the  "tunneling  direction",  at  which  time  a  WKB  tunneling 
probability  is  computed.  Makri  and  Miller^"^  assume  a  more  general  tunneling  coordinate, 
specificalh’.  the  shortest  straight-line  distance  connecting  the  caustics  of  the  (assumed 
quasiperiodic  and  "quantized")  trajectories  in  the  reactant  and  product  wells.  Thus,  in  general,  the 
Makri-Miller  tunneling  direction  will  correspond  to  motion  of  all  the  atoms  in  the  system.  In  this 
system,  a  tunneling  coordinate  defined  in  terms  of  a  single  internal  coordinate  (a  la  Waite  and 
Miller^^)  is  not  adequate;  passage  from  the  equilibrium  geometry  to  the  transition  state  requires 
significant  displacement  of  both  the  bond  and  angle  coordinates.  Thus,  the  more  general  tunneling 
direction  described  by  Makri  and  Miller  is  indicated. 

Rigorously,  Makri  and  Miller^-^  require  that  the  classical  trajectories  be  quasiperiodic  and 
correspond  to  a  semiclassical  state  of  the  system.  The  latter  can  be  accomplished  by  using 
adiabatic  switching.^5  We  did  not  attempt  a  semiclassical  quantization  of  the  system  but,  rather, 
chose  to  investigate  the  influence  of  various  ad  hoc  distributions  of  energy  on  the  semiclassical 
tunneling  splitting.  Initial  conditions  for  the  trajectories  were  obtained  by  partitioning  the  ground- 
state  energy  using  four  different  schemes.  Specifically,  the  energy  was  partitioned  i)  entirely  into 
the  zeroth-order  stretch  or  ii)  the  bending  mode;  iii)  in  accordance  with  standard  quasiclassical 
prescriptions;  or  iv)  by  totally  random  assignment  of  energy  (with  respect  to  the  Canesian 
coordinates  and  momenta). 


The  calculated  tunneling  splittings  are  given  in  Table  1.  There  is  a  clear  dependence 


Table  1.  Calculated  tunneling  splittings 

Excitation  Scheme 

Calculated  Splitting  (cm'*) 

Pure  angle 

16.2 

Pure  stretch 

41.5 

Random 

38.8 

Quasiclassical 

50.9 

on  the  nature  of  the  initial  conditions,  with  the  computed  splitrings  ranging  over  roughly  a  factor  of 
three.  Interestingly,  the  results  compare  quite  favorably  to  the  measured  value^^  of  2 1  cm'  1 .  The 
largest  splitting,  50.9  cm'^.  corresponds  to  quasiclassical  selection.  Tlie  smallest  splitting.  16,2 
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cm'^  is  obtained  for  pure  angle-bending  excitation.  Pure  bond-stretching  and  totally  random 
excitation  yield  splittings  of  41.5  and  38.3  cm'^,  respectively. 

The  results  shown  in  Table  1  are  qualitatively  reasonable.  The  zero-point  energy  (6.6 
kcal/mol)  classical  dynamics  are  essentially  quasiperiodic.^'  Projection  of  all  of  the  energy  onto 
the  bending  coordinate  results  in  a  stable,  highly  excited  bending  motion  and  very  little  energy 
tends  to  leak  into  the  stretching  mode.  The  tunneling  points  are  "far"  from  the  transition  state.  By 
contrast,  localized  excitation  of  the  stretching  mode  leads  to  a  set  of  turning  points  that  are  closer  to 
equilibrium,  which  means  that  the  "distance"  from  the  reactant  well  to  the  product  well  is 
decreased;  and  the  height  of  the  barrier  through  which  tunneling  must  occur  is  decreased  as  well. 
These  two  factors  lead  to  a  much  increased  splitting  for  the  case  of  pure  stretching  excitation. 
Quasiclassical  excitation  partitions  energy  into  both  the  bending  and  stretching  modes,  suggesting 
that  the  tunneling  splitting  will  be  larger  than  for  either  of  the  two  more  selective  excitations. 

Totally  random  distribution  of  the  energy  spans  a  considerably  larger  volume  of  phase  space  than 
the  other  three  methods.  Thus,  it  is  likely  that  some  trajectories  have  very  large  probabilities  of 
tunneling  (leading  to  a  large  splitting)  while  others  will  have  correspondingly  small  tunneling 
probabilities. 

The  intuitive  arguments  of  the  preceding  paragraph  are  illustrated  by  tuialysis  of  plots  of  the 
the  tunneling  points  of  the  tunneling  paths  for  ensembles  of  trajectories.  The  reason  for  the 
reduced  splitting  for  pure  angle-bending  excitation  is  that  the  turning  points  for  tunneling 
trajectories  are  tightly  clustered  such  that  the  distance  to  the  product  "state"  in  the  product  well  is  on 
the  order  of  0.8  bohr,  and  the  height  of  the  barrier  which  must  be  penetrated  is  significant.  Thus, 
tunneling  must  occur  through  a  tall,  wide  potential  -  and  the  probability  of  tunneling  falls  of 
exponentially  with  the  area  of  the  barrier  which  must  be  penetrated.  The  situation  for  pure 
stretching  excitation  is  more  conducive  to  tunneling.  First,  the  average  potential  energy  at  the 
tunneling  points  is  higher.  Consequently,  the  average  height  of  the  effective  barrier  is  decreased. 
Also,  the  width  of  the  barrier  is  smaller.  A  final  factor  is  the  frequency  of  the  motion  in  the 
tunneling  direction.  On  average,  four  tunneling  points  were  encountered  for  pure  bending 
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trajectories.  By  contrast,  there  were  10-11  per  pure  bond-stretching  trajectory.  Thus,  not  only  is 
the  probability  of  penetrating  the  barrier  at  each  turning  point  greater  for  the  case  of  bond-stretching 
excitation,  but  the  frequency  of  turning  points  is  more  than  twice  as  large.  Quasiclassical  excitation 
results  in  a  broad  swath  of  turning  points.  The  individual  tunneling  probabilities  fall  off  rapidly 
with  increasing  values  of  the  ordinate.  Individual  quasiclassical  trajectories  typically  exhibit  a  wide 
range  of  individual  tunneling  probabilities  due  to  the  relatively  equitable  partitioning  of  energy 
between  the  two  modes.  Totally  random  sampling  allows  for  exploration  of  a  much  larger  volume 
of  the  potential.  The  result  is  that  the  typical  trajectory  is  less  likely  to  explore  regions 
corresponding  to  high  tunneling  probabilities. 

The  results  suggest  that  it  may  be  possible  to  observe  mode-specific  effects  in  the  tunneling 
splitting  of  real  molecules.  Since  tunneling  is  a  sensitive  function  of  the  dynamics  and  the  energy, 
selective,  low-energy  excitations  could  perhaps  be  designed  to  excite  selected  modes  (i.e.,  C-O-H 
bends  or  skeletal  modes)  in  the  hopes  of  observing  mode-specific  splitting.  Since  the  principal 
interest  would  be  on  low  energies,  the  description  could  probably  be  couched  in  terms  of  normal 
modes  and,  moreover,  the  lifetime  of  the  initially  prepared  states  might  be  long  enough  to  be 
chemically  significant. 

In  summary,  this  study  of  the  tunneling  splitting  for  a  model  of  the  malonaldehyde 
molecule  indicates  significant  mode-specific  effects.  The  mode  specificity  is  understandable  in 
terms  of  the  quasiperiodic  dynamics  that  occur  on  the  potential-energy  surface.  It  is  imponant  that 
practical  methods  be  developed  for  treating  tunneling  effects  in  large  systems,  and  the  semiclassical 
approach  we  have  used  offers  the  possibility  of  a  general  method. 

The  simplest  case  of  tunneling  is  where  the  "tunneling  coordinate"  can  be  defined  (at  least, 
approximately)  as  a  normal  mode.  We  have  carried  out  a  semiclassical  calculation^^  of  the  effects 
of  vibrational  energy  on  tunneling  in  the  cis-trans  isomerization  in  HSiOH  using  a  potential-energy 
surface  based  on  the  available  experimental  and  theoretical  data.  The  isomerization  channels.  HOSi 
bend  and  torsion,  both  have  barriers  on  the  order  of  1 1  kcal/mol.  The  effect  of  internal  rotation 
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increases  the  barrier  along  the  bending  mode.  The  calculations  show  strong  mode-specific 
behavior  in  the  tunneling  rate. 

The  calculated  tunneling  rates  as  a  function  of  the  excitation  energy  show  that  when 
the  total  energy  of  the  system  is  low  (15  kcal/moi  or  lower),  the  computed  tunneling  rate 
coefficient  is  three  orders  of  magnitude  larger  for  the  cis-^trans  than  that  for  the  trans 
—>cis.  The  ratio  of  the  tunneling-corrected  rate  to  the  purely  classical  rate  increases  with 
decreasing  total  energy. 

C.  Mode  Selective  Dynamics  in  Highly  Vibrationally  Excited  Polyatomic 

Molecules 

Numerous  studies  have  shown  that  the  dynamics  of  few-mode  models  and  molecules  are 
chaotic  at  "high  energies."  33a,36-54  However,  it  is  not  as  clear  to  what  extent  chaos  is  a  factor  in 
unimolecular  reactions  involving  large  molecules.54.55  Chaos  obviously  leads  to  statistical 
behavior  in  chemical  reactions.  The  validity  of  statistical  theories  (such  as  RRKM)  and  the 
possibility  of  observing  mode- selective  chemistry  depend  on  the  nature  of  the  dynamics  at  energies 
at  and  just  above  reaction  thresholds.  While  a  large  body  of  literature  deals  with  the  phase-space 
structure  of  small  (mostly  three  atoms),  highly  excited  molecules,56  relatively  little  has  been 
reponed  for  the  dynamics  of  large  molecules.  This  is  primarily  due  to  difficulties  in  extending  the 
theoretical  techniques  used  to  study  small  molecules  to  large-molecule  dynamics. 

Among  the  more  common  ways  of  "detecting"  chaos  are  techniques  such  as  monitoring  the 
separation  in  phase  space  of  "initially  adjacent"  trajectories40'43,44  or  examination  of  Poincaire 
surfaces  of  section.44 .46,57  Brumer  and  Duff40  have  shown  that  the  rate  of  divergence  of  initially 
adjacent  trajectories  can  be  used  as  a  tool  for  diagnosing  chaotic  dynamics.  Carter  and  Brumer42 
and  Davis  and  Wagnei45,44  published  a  number  of  papers  using  this  method  to  study  chaotic 
vibrational  motion  in  OCS.  The  surface-of-section  method  is  usually  restricted  to  two  degrees-of- 
freedom  systems  although  generalizations  to  3-D  systems  have  been  developed.^^a  Another 
simple,  though  not  totally  satisfactoi.  way  of  studying  the  dynamics  of  large  systems  is  by 
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calculating  power  spectra  from  classical  trajectories.37.42.45,46.54b,58-65  Power  spectra  can 
provide  information  about  the  structure  of  the  phase  space.  For  instance,  purely  periodic  motion 
gives  a  "stick"  spectrum.  By  contrast,  chaotic  motion  gives  a  broadband  spectrum  having  a 
characteristic  pattern  that,  strictly  speaking,  only  emerges  in  the  long-time  limit.^^’*^^  From  a  more 
pragmatic  point  of  view,  a"grassy"  spectrum  is  usually  associated  with  chaos  in  the  (relatively) 
shon-time  dynamics  of  small  systems.  For  most  conditions  of  interest  the  spectra  will  likely  be 
intermediate  between  these  two  extremes,  and  will  be  such  that  precise  conclusions  cannot  be 
drawn.  Nevertheless,  power  spectra  are  useful  for  obtaining  information  about  the  qualitative 
behavior  of  the  dynamics. 

Among  the  earliest  publications  in  the  chemical  dynamics  literature  concerning  power 
spectral  analysis  of  the  vibrational  dynamics  of  "chemical"  systems  is  a  paper  by  Noid, 
Kozykowski,  and  Marcus.^^  Since  then,  there  have  been  numerous  studies  that  used  power 
spectra  as  a  qualitative  means  of  classifying  the  classical  dynamics  as  being  either  regular  or 
chaotic.^^-'^^’^^’^^'^ 

Smith  and  Shirts^^  concluded  on  the  basis  of  power  spectra  that  the  dynamics  of  HNC 
isomerization  become  chaotic  at  a  fairly  well-defined  energy.  They  reported  "grassy "  structure  and 
broadening  in  the  power  spectra.  Averaging  the  spectra  over  a  small  ensemble  of  trajectories 
tended  to  smooth  out  the  grassiness  but  the  broad  features  remained  prominent  in  the  spectra. 
Sumpter  and  Thompson^^-^^  used  power  spectra  to  qualitatively  examine  the  effects  of  rotational 
motion  on  the  vibrational  dynamics  of  H2O2.  The  presence  of  angular  momentum  results  in  a 
substantial  broadening  in  the  peaks  relative  to  the  nonrotating  case  and  introduces  new  peaks  into 
the  spectrum. 

We  have  employed  power  spectra  in  an  attempt  to  ascertain  whether  the  dynamics  of 
several  polyatomic  molecules  are  primarily  chaotic  or  nonchaotic.^^  The  molecules  we  chose  to 
study  are:  SiFa,  C2H4,  CH3ONO,  and  •CH2CH2CI.  The  results  indicate  the  existence  of  well 
resolved  (and  assignable)  bands  in  the  vibrational  spectrum  of  SiFa,  CH3ONO,  and  (to  a  lesser 
extent)  •CH2CH2CI  for  energies  up  to  and  above  the  threshold  for  unimolecular  dissociation. 
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These  bands  are  found  to  be  present  at  low  energies  (i.e.,  the  zero-point  energy)  and  remain 
relatively  unchanged  over  a  wide  range  of  internal  energies.  By  contrast,  the  vibrational  spectrum 
of  C2H4  indicates  a  transition  to  chaotic  motion  at  energies  below  the  reaction  threshold.  The 
implication  is  that,  contrary  to  conclusions  often  drawn  for  "small"  molecular  systems,  the 
dissociation  of  large  polyatomic  molecules  is  likely  not  to  be  governed  by  chaotic  dynamics. 

We  have  used  classical  trajectories  on  an  ab  i/iino-based  potential-energy  surface  to  study 
the  intramolecular  dynamics  and  unimolecular  dissociation  of  the  CF2Br  radical.^^  Our  interest 
here  was  to  study  mode  localization  and  its  influence  on  reaction.  The  CF2Br  radical  has  been 
observed  as  a  primary  product  of  the  248  nm  photolysis  of  CF2Br2.^®’^^  Krajnovich,  Zhang. 
Butler,  and  Lee^^  examined  the  photochemistry  of  CF2Br2  in  crossed  laser-molecular  beams. 

They  concluded  that  the  only  primary  reaction  channel  is  the  elimination  of  a  Br  atom.  They  also 
detected  CF2  radicals  which  they  attributed  to  efficient  secondary  photodissociation  of  the  CF^Br 
radical.  Gosnell,  Taylor,  and  Lyman^^  used  ultrafast  time-resolved  absorption  spectroscopy  to 
study  the  photodissociation  of  CF2Br2.  They  also  detected  the  CF2Br  radical  as  a  primary 
product.  However,  they  concluded  that  the  CF2  radical  is  formed  in  the  ground  or  first-excited 
vibrational  state  of  its  V2  bending  mode  during  the  "unimolecular  reaction  of  a  vibrationally  hot 
CF2Br  intermediate  species".  In  addition,  they  concluded  that  the  evolution  of  the  CF2Br 
absorption  spectrum  over  6  ps  is  indicative  of  slow  intramolecular  vibrational  redistribution.  These 
results  are  sufficiently  intriguing  to  warrent  a  study  of  the  dynamics  of  this  system. 

Power  spectra  of  CF2Br  show  strong  evidence  of  mode  localization  .  Noid, 

Koszykowski,  and  Marcus^^  examined  power  spectra  for  two  model  Hamiltonians  and  attributed 
sharp  line  spectra,  for  which  peaks  are  assignable  as  fundamentals,  overtones,  and  combinations, 
to  quasiperiodic  dynamics.  In  contrast,  they  demonstrated  that  when  the  dynamics  are  known  to 
be  ergodic  that  the  spectra  become  grassy  and  broadened.  However,  as  noted  by  Noid, 
Koszykowski,  and  Marcus,  centralization  of  certain  frequencies  is  still  apparent  in  the  "chaotic 
spectra"  indicative  of  some  form  of  mode  localization  in  the  frequency  domain. 
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Tlie  sharp  peaks  in  a  regular  power  spectrum  will  persist  with  increasing  time  with 
intensities  which  increase  with  time  and  widths  which  are  inversely  proportional  to  time  (due  to  the 
increase  in  spectral  resolution  with  increasing  time),  whereas  chaotic  power  spectra  become  more 
grassy  and  do  not  display  peaks  which  persist  with  increasing  time.^^  This  is  well  demonstrated  in 
a  study  by  Smith  and  Shins^^  of  HCN.  They  showed  that  the  ensemble  averaged  spectra  do  not 
exhibit  the  expected  grassiness  with  increasing  energy  but  that  individual  power  spectra  do. 
However,  the  peaks  in  the  ensemble  averaged  spectra  decreased  in  intensity  and  broadened,  as 
expected,  as  the  number  of  bending  mode  quanta  was  increased.  Dumont  and  Brumer^-'^  carried 
out  a  detailed  analysis  of  the  relationship  between  sharp/grassy  power  spectra  and  regular/chaotic 
dynamics.  They  concluded  that  "individual  power  spectra,  computed  over  a  finite  time,  are  unable 
in  principle  to  provide  a  distinction  between  formally  regular  and  formally  chaotic  dynamics"  but 
that  "if  one  adopts  the  view  that  formal  definitions  of  chaotic  vs  regular  dynamics  are  less  relevant 
than  behavior  over  a  finite  time  scale  of  observational  interest,  then  simple  vs  complex  power 
spectra  do  provide,  by  definition,  useful  distinguishing  characteristics".  It  is  the  second  view 
which  we  have  adopted  in  our  analysis  of  the  power  spectra  which,  when  examined  in  conjunction 
with  autocorrelation  functions  of  the  internal  coordinates,  provide  evidence  for  mode  localization  in 
CF2Br. 

Five  excitation  energies,  measured  relative  to  the  zero-point  energy  of  the  harmonic  normal 
modes  (approximately  6  kcal/mol),  w-ere  considered.  The  lowest  excitation  energy  is  35  kcal/mol 
which  is  slightly  higher  than  the  C-Br  bond  well-depth  of  33.3  kcal/mol.  Krajnovich  et 
estimated  the  maximum  product  internal  energy  of  CF2Br  to  be  39  kcal/mol.  Other  excitation 
energies  (45,  50,  and  60  kcal/mol)  were  considered  to  determine  what  effect  (if  any)  higher 
excitation  energies  might  have  on  the  rate.  Seven  calculations  were  carried  out  at  each  of  the  five 
energies.  Six  of  these  calculations  involved  mode-specific  excitation  of  the  individual  normal 
modes.  Tlie  mode-specific  initial  conditions  were  determined  by  populating  all  but  one  of  the 
normal  modes  with  its  calculated  ground-state  vibrational  energy.  The  remaining  normal  mode 
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was  assigned  a  total  energy  which  is  the  sum  of  the  zero-point  energy  for  that  mode  plus  an 
arbitrary  excitation  energy . 

Ab  iniiio  calculations  were  performed  on  CF2Br  at  the  Hartree-Fock  level  using  the  MESA 
electronic  structure  codes^^  with  6-3 IG  basis  sets^^  for  C  and  F  and  a  DZ  basis  on  Br  where  an 
effective  core  potential^'^  was  employed  to  replace  the  (Arl3d*0  core.  The  geometry  was  optimized 
using  analytical  gradient  techniques  and  the  vibrational  frequencies  were  determined  at  the  optimum 
geometry  using  numerical  second  derivatives.  The  optimized  geometrical  parameters  were  as 
follows:  R(C-Br)=1.953  A,  R(C-F)=1.346  A,  0(Br-C-F)=l  14.3‘.  and  e(F-C-F)=l  10.4'.  For 
comparison,  the  calculated  C-F  bond  length  in  CF2  is  1.335  A  in  the  6-3  IF  basis  used  here, 

1 .283  A  in  the  6-3 IG*  basis,  and  1.304  A  experimentally.^^  This  information,  along  with 
experimental  data,  was  used  to  formulate  a  global  potential-energy  surface. 

The  computed  power  spectra  for  CF2Br  clearly  demonstrate  mode  localization.  There  is  a 
sharp  peak  at  about  the  vg  (Br-C-F  bend)  frequency  in  all  the  power  spectra  for  the  various  initial 
conditions.  Furthermore,  in  a  series  of  power  spectra  for  5  ps  time  windows  of  a  45  ps  time 
history,  the  peak  remains  well-defined  over  the  entire  time  history.  The  power  spectra  show  that  at 
excitation  energies  well  above  the  C-Br  bond  dissociation  energy  the  dynamics  are  not  globally 
chaotic  and  that,  in  fact,  the  frequency  bands  due  to  particular  normal  coordinates  are  still 
identifiable.  However,  the  calculated  rate  constants  show  that  sufficient  interaction  must  be 
occurring  between  vg  and  other  modes  to  prevent  a  decrease  in  the  rate  of  reaction  when  initially 
contains  all  of  the  excitation  energy.  The  rate  constants  demonstrate  mode  selectivity  only  with 
respect  to  direct  excitation  of  the  reaction  coordinate  mode  (V4);  there  is  no  selectivity  in  the  rates 
for  excitations  of  the  other  normal  modes.  Furthermore,  except  for  V4,  the  rate  constants  are 
comparable  to  those  for  equipartition  of  the  excitation  energy  among  the  normal  modes. 

These  results  demonstrate  that  for  the  CF2Br  Hamiltonian  used  in  this  study,  mode 
localization  is  not  a  sufficient  condition  for  mode-specific  reaction  rates.  Additional  studies 
employing  state-specific  initial  conditions  would  be  useful  for  examining  the  effect  of  mode 
localization  on  the  rates  of  chemical  reactions. 
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As  illustrated  by  the  above  study  and  many  others,  mode  selectivity  is  not  genenilK’ 
observ'ed  at  the  energies  required  to  break  a  chemical  bond.  In  contrast,  the  vibrational 
predissociation  of  van  der  Waals  systems  exhibits  non-statistical  behavior.  Although  van  der 
Waals  molecules  are  distinctly  different  from  chemically  bound  systems,  they  serve  as  useful 
examples  of  mode  selective  systems  and  help  to  illustrate  some  fundamental  aspects  of 
nonstatistical  behavior.  In  an  effort  to  understand  mode  selectivity  in  general,  it  is  useful  to  study 
systems  that  are  intermediate  between  van  der  Waals  and  chemically-bound  molecules.  In  this 
regard  we  have  reponed  a  series  of  intramolecular  conversions  that  have  energy  barriers  on  the 
order  of  10-15  kcal/mol.^^  More  recently,  we  studied  a  "van  der  Waals/chemically  bound" 
molecule,  N2O3. 

We  used  classical  trajectories  to  study  mode  selective  vibrational  predissociaiion  in 
N203.^7  The  results  are  in  accord  with  the  conclusion  of  Chewter,  Smith,  and  Yarwood^^  that 
N2O3  vibrational  predissociation  exhibits  non-RRKM  behavior.  In  fact,  we  find  strong  mode 
selective  behavior  as  evidenced  by  the  variation  in  the  rate  as  a  function  of  the  initial  panitioning  of 
energy  among  the  normal  modes.  This  variation  in  the  rates  can  be  explained  by  the  extent  to 
which  the  excited  modes  project  onto  the  reaction  coordinate.  The  nitro  group  bending  (V4)  and 
symmetric  stretching  (V3)  modes  involve  motion  with  significant  components  directed  along  the  N- 
N  bond  (approximately  the  reaction  coordinate),  while  the  nitroso  group  stretch  (vj)  and  the  niiro 
group  asymmetric  stretch  (V2)  do  not.  The  trajectory  results  show  that  V4  mode  excitation  is  the 
most  efficient  in  promoting  dissociation  and  V3  is  the  second  most  effective  mode.  The  lifetimes 
for  V]  and  V2  excitations  are  significantly  longer  than  those  for  V3  and  V4  excitations. 

D,  The  "Zero-Point  Energy  Problem"  in  Classical  Trajectories 

Classical  trajectories  have  been,  and  continue  to  be.  a  crucial  component  of  the  "toolkit”  in 
chemical  dynamics.  Classical  simulations  are  feasible  for  many  problems  where  quantum 
treatments  are  not.  Fonunately,  many  chemical  processes  are  essentially  classical.  However, 
some  systems  are  quantum  mechanical  and  there  are  situations  where  classical  mechanics  fail.  One 
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such  problem  is  the  failure  of  classical  mechanics  to  preserve  zero-point  energy  constraints.  We 
have  raised  the  question/''^  of  whether  it  is  possible  to  "fix"  the  zero-point  energy  problem  w  iihout 
altering  the  dynamics  such  that  equally  or  more  serious  aphysical  behavior  is  introduced.  Various 
approaches  have  been  tried  to  correct  this  problem.^0-85 

Some  studies  have  been  done  to  determine  the  contribution  of  the  zero-point  energy  to 
reaction  rates  in  classical  studies  by  including  only  a  fraction  of  the  zero-point  energy  of  the 
molecule.  Guan  et  al.^^  used  this  approach  in  studies  of  the  near-threshold  ovenone-excitation 
induced  photodissociation  of  H2O2.  They  found  that  there  is  a  significant  dependence  of  the 
classically  calculated  dissociation  rates  on  the  amount  of  zero-point  energy.  However,  they  found 
that  mode-specific  effects  are  qualitatively  the  same  for  various  amount  of  zero-point  energy.  Lu 
and  HaseSh82  studied  the  effect  of  zero-point  energy  in  calculations  of  the  linewidth  of  C-H 
ovenones  following  excitation  of  a  C-H  local  mode  in  benzene.  They  added  only  a  fraction  of  the 
zero-point  energy  (0  ^  fzpc  ^  1)  to  the  modes  of  the  system  and  they  excluded  zero- point  energy 
from  selected  normal  modes.  They  found  that  excluding  zero-point  energy  from  7  normal  modes 
having  strong  CCH  bend  character  had  a  significant  effect  on  the  linewidth.  By  contrast, 
excluding  zero-point  energy  from  the  other  14  normal  modes  but  adding  it  to  the  7  CCH  modes 
gave  results  in  near  agreement  with  the  full-ZPE  case,  except  for  higher  overtone  excitations.  Lu 
and  Hase  also  considered  the  effects  of  decreasing  the  number  of  degrees  of  freedom  in  the  system 
in  an  effon  to  reduce  the  effects  of  aphysical  classical  energy  flow.  In  particular,  they  extracted 
ponions  of  the  force  field  for  the  planar  benzene  molecule  (C6H6)  to  yield  force  fields  for  two 
fragments  (C3H  and  C3H3)  and  compared  linewidths  for  the  two  fragments  with  experimental 
measurements.  The  effect  of  doing  so  was  to  decrease  the  total  energy  content  of  the  system  and 
change  the  intramolecular  couplings.  They  found  that  the  overtone  linewidths  in  the  fragments 
were  in  considerably  better  agreement  with  the  experimental  measurements  for  benzene  excited  to 
the  corresponding  state. 

Another  possible  route  to  understanding  the  effects  of  zero-point  energy  flow  in  classical 
trajectory  calculations  is  illustrated  by  the  work  of  Nyman  and  Davidsson^^  on  the  exchange 
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reaction  0{-'’P)  +  OHi^FI  )  02(^1’';,')  +  H(2S).  They  staned  with  "all-or-none  "  of  the  zero-point 

energy  in  the  OH(2/7)  vibration,  and  then  selectively  "threw  out"  trajectories  which  did  not  saiisfv 
predefined  criteria  regarding  zero-point  energy  in  the  reaction  products  (02(^2''g)  or,  in  the  case  of 
unreactive  scattering,  OH(2/7)].  The  critera  range  from  retaining  all  trajectories  regardless  of  the 
final  vibrational  energy  to  retaining  only  those  trajectories  that  had  at  least  the  zero-point  energy. 
The  results  showed  a  strong  dependence  on  which  scheme  of  categorizing  the  trajectories  is 
employed;  The  treatment  that  yields  the  best  agreement  between  the  calculated  and  experimental 
results  depends  upon  what  properties  are  being  compared. 

These  examples  serve  to  exemplify  "passive"  methods  for  dealing  with  the  zero-point 
energy  problem  and  also  illustrate  the  effects  of  zero-point  energy  flow.  They  are  based  on  the  idea 
of  allowing  uie  trajectories  to  evolve  according  to  the  classical  equations  of  motion.  By  contrast, 
an  "active"  method  for  controlling  the  flow  of  zero-point  energy  during  the  course  of  individual 
trajectories  has  been  proposed  simultaneously  by  Bowman,  Gazdy,  and  Sun^^  and  Miller,  Hase. 
and  Darling.^5  The  approach  by  Bowman  er  al.  is  encompassed  in  that  of  Miller  et  al.  The  gist  of 
both  methods  is  that,  so  long  as  none  of  the  (assumed  separable)  mode  energies  in  a  system  fall 
below  the  zero-point  level,  it  should  be  allowed  to  evolve  according  to  the  classical  equations  of 
motion.  However,  when  one  or  more  of  the  mode  energies  falls  below  the  zero-point,  an 
instantaneous,  external  force  is  applied  so  as  to  exactly  reverse  the  velocity  in  the  models)  in  \\  hich 
the  energy  has  dropped  below  the  zero-point  level.  If  the  kinetic  energy  is  diagonal,  then  the  total 
energy  will  be  conserved.  (However,  the  angular  momentum  will  not  be  conserved  except  in 
special  cases  such  as  linear  systems.) 

We  have  addressed  the  question  of  whether  or  not  it  is  feasible  to  employ  "active" 
constraints  to  deal  with  the  zero-point  energy  problem.^^  To  this  end  we  apply  the  methods 
suggested  by  Bowman  ei  and  Miller  et  to  the  simplest  situation  one  can  imagine,  the 
dynamics  of  the  ground  state  of  a  system  of  coupled  oscillators  (specifically,  the  Henon-Heiles 
Hamiltonian).  The  results  given  by  the  usual,  i.e.,  unconstrained  classical  mechanics  and  by  the 
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zero-point  energy  constrained  classical  mechanics  are  compared  to  the  correct  quantum  mechanical 
calculation  of  several  quantities  of  interest. 

We  find  that  for  this  system  of  coupled  oscillators  (Henon-Heiles)  initially  in  its 
(approximate)  ground  state,  there  is  no  zero-point  energy  problem  as  illustrated  by  a  comparison  of 
quantum  and  classical  results.  We  find  that  the  proposed  model^'^'^5  for  correcting  problems  with 
zero- point  energy  gives  mode  energies  somewhat  is  better  agreement  with  the  quantum  results,  but 
the  unconstrained  classical  trajectories  are  not  in  serious  error  either.  This  difficulty  arises  when 
one  examines  power  spectra  and  plots  of  the  configuration  space  projections.  The  unconstrained 
classical  mechanics  gives  a  single  peak  in  the  power  spectra,  indicating  quasiperiodic  motion, 
while  the  constrained  classical  mechanics  spectra  show  broad,  essentially  continuous  distributions 
of  frequencies,  characteristic  of  chaotic  motion.  The  configuration  space  plots  confirm  the  chaotic 
behavior  in  the  constrained  trajectories.  The  quantum  results  (spectra)  are  clearly  closer  to  those  of 
the  unconstrained  classical  mechanics.  We  examined  a  variety  of  initial  conditions,  but  all  the 
qualitative  results  are  as  described  above.  The  constrained  dynamics  very  effectively  prevent  the 
energy  in  each  mode  from  falling  below  its  local  zero  point  value,  but  does  severe  "damage "  to  the 
dynamics.  In  other  words,  the  constraints  severly  modify  the  classical  dynamics  in  a  physically 
incorrect  fashion.  Thus,  our  work  (nor  that  of  others  as  far  as  we  know)  has  not  provided  general 
solution  for  this  vexing  problem.  We  remain  interested  and  plan  to  pursue  it.  Perhaps  the 
problem  with  the  above  approach  is  that  the  "hard  wall"  character  of  the  constraint  is  too  harsh, 
inducing  chaotic  behavior  in  the  dynamics  and  a  "kinder,  gentler"  version  of  the  model  can  be 
devised  (e.g.,  a  soft  repulsive  interaction  in  the  action  variable),  but  it  remains  to  be  seen  whether 
any  strictly  classical  "fix"  can  be  successful  and  useful. 

One  may,  of  course,  decide  to  give  up  on  a  strictly  classical  approach  and  try  to  utilize 
various  semiclassical  approximations.  Alimi  et  al.^  have  recently  suggested  an  approximate 
version  of  Heller's^^  semiclassical  wavepacket  methods  for  this  purpose.  This  is  a  promising 
approach  and  certainly  wonh  pursuing,  though  one  notes  that  it  involves  more  effort  to  apply  than 
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does  a  purely  classical  treatment  (e.g.,  matix  elements  of  the  potential  energy  surface  must  be 
computed),  and  Gaussian  wavepackets  can  give  unphysical  results  in  very  anharmonic  systems. 
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